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Abstract: The single-species reaction-diffusion process A+A — > O is examined in the presence 
of an uncorrelated, quenched random velocity field. Utilising a field-theoretic approach, we find 
that in two dimensions and below the density decay is altered from the case of purely diffusing 
reactants. In two-dimensions the density amplitude is reduced in the presence of weak disorder, 
yielding the interesting result that Sinai disorder can cause reactions to occur at an increased 
rate. This is in contrast to the case of long-range correlated disorder, where it was shown that 
the reaction becomes sub-diffusion limited. However, when written in terms of the microscopic 
diffusion constant it is seen that increasing the disorder has the effect of reducing the rate of the 
reaction. Below two dimensions, the effect of Sinai disorder is much more severe and the reaction 
is shown to become sub-diffusion limited. Although there is no universal amplitude for the time- 
dependence of the density, it is universal when expressed in terms of the disorder-averaged diffusion 
length. The appropriate amplitude is calculated to one-loop order. 
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1 Introduction 



The dynamical, many-body problem of diffusing and reacting chemicals provides an ideal testing 
ground for the many methods of non-equilibrium physics. Such reaction systems occur in nature 
in a wide variety of guises, from conventional chemical reactions to more exotic processes such as 
domain coarsening in magnetic systems and exciton annihilation in crystals ||]. When an 

equilibrium has been reached in a reversible reaction, the methods of statistical mechanics can 
be used to give the appropriate reactant-density ratios. However, in the decay to equilibrium Q 
and in irreversible processes no such unified approach yet exists. Nevertheless, over the last few 
decades a great many studies have been made utilising numerical, exact and renormalisation-group 
techniques, and the time-dependence of reaction systems has been found to be remarkably rich. 
It has been shown that the collective behaviour of the reaction system is rather sensitive to the 
statistical properties of a single reactant's motion, with a basic example given by the case of the 
single-species diffusive reaction A + A — > O. This is because the re-entrancy of random walks in 
two dimensions and below alters the density decay from the simple mean-field prediction. In these 
lower dimensions, the density loses all dependence on the microscopic reaction rate and becomes 
a universal function of only the diffusion length: the reaction is diffusion limited. In this paper 
we examine the behaviour of the A + A — > O reaction process in the presence of an uncorrelated, 
quenched random velocity field. Before going on to describe our results and method, we review 
two cases, first the pure reaction-diffusion process, and second some studies of the reaction process 
in correlated long-range potential disorder. 

The mutually annihilating random walk (MARW) A + A — > O is a fundamental theoretical 
model in the study of reaction systems. It describes the process whereby diffusing A particles 
may react pairwise, at a rate A on contact. The MARW and coalescing random walk A + A — > A 
are members of the same universality class |J, and therefore the results below are also valid for 
the coalescence process (albeit with trivial changes in the prefactors). The basic starting point 
in the analysis of such systems is the mean-field or rate equation. This corresponds to writing a 
self-consistent equation for the average reactant density n as a function of t, that ignores all spatial 
correlations (and more seriously anti- correlations) , 

drx 1 

— - = DV 2 n — 2Xn 2 , with the late time result n ~ — — (1) 

dt 2\t y ' 

where we have introduced the diffusion constant D and the reaction rate is A. Neglecting the 
effect of correlations is equivalent to assuming that the reactants remain well mixed throughout 
the reaction process. However, due to the statistics of random walks in two dimensions and below, 
simple diffusion of particles itself is not sufficiently fast to maintain a well-mixed state. For this 
reason, the mean-field result (|l|) loses its validity in two dimensions and below. In these lower 
spatial dimensions the re-entrancy of random walks means that reactants come into contact many 
times, each time providing an opportunity for a reaction to take place. This implies that even a 
small reaction rate A does not limit the global rate of reaction. Therefore, in two dimensions and 
below the reaction becomes diffusion limited with the density decays 

n= logCDt) d = 2 n ~_^!_ ford<2 (2) 

8irDt (Dt) d / 2 v ; 

where Ad is a universal amplitude, a function only of the dimension d. This amplitude has been 
calculated via an e expansion for d = 2 — e dimensions || (also given in equation (|8|)) and the 
exact result Ai = (87r) -1 / 2 for one dimension can be found in reference Q. 
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Over the last few years there has been an increased interest in the behaviour of reaction systems 
with reactants that perform motion different from pure diffusion. As well as ballistic gas-phase 
reactions [||-[1C], studies have appeared with reactants that perform diffusion in the presence of 
turbulence [11 1 12] and also in quenched random velocity fields. Some basic categories for the 
statistical properties of these random velocity fields have been identified, see for example jnj [14]. 
In particular, a distinction can be made between uncorrelated Sinai disorder and the long-range 
correlated potential disorder, with the momentum-space correlator 7/fc 2 . Studies have been made 
of the behaviour of the reaction front in the segregated two-species A+B — ► O reaction with various 
forms of Sinai disorder [|D|]-[17] and correlated potential disorder [18] 19]. A comprehensive study 
has been made of the A + A — > O scheme with correlated potential disorder (2(| 21] in two 
spatial dimensions and solutions also exist for this single-species reaction with random barriers 
and random traps [p2]] . Furthermore, the specific case of single-species reactions in Sinai disorder 
was recently examined in one dimension in the context of aging phenomena, with the persistence 
exponent derived [23|. In the two-dimensional case of A + A — > O in potential disorder, it was 
found that the reaction process becomes sub-diffusion limited: the reaction rate is lower than the 
MARW. The form for the density with this potential disorder, which should be compared with the 
two-dimensional MARW in equation (0), was found to be 



n 



X*t 
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with 5 
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> 0, A* = 3£>/? 2 7 



(3) 



where f3 2 j measures the disorder strength, 5 is a non-universal exponent, and A* is an effective 
reaction constant. The interpretation was that the long-range disorder produces potential traps 
on all length scales which, after some time, will contain at most one reactant. Reactions can then 
only occur when the trapped reactants move between traps. However, at the same time, as the 
reactants explore the landscape, they get caught by increasingly deep traps, leading to sub-diffusive 
motion. 

In this paper, we present results from a study of the single-species A + A — > O reaction process 
in the presence of an uncorrelated, quenched random velocity field in dimension two and below. 
Despite the lack of long-range correlations the kinetic behaviour of diffusing particles is changed, 
as can be seen in the behaviour of the diffusion-length squared (r 2 ). In two dimensions, it is 
known from a renormalisation-group (RG) treatment that this length is altered by the presence of 
a logarithm [24]. Below two dimensions, d = 2 — e, the RG gives the two-loop result for the dynamic 
exponent of z = 2 + 2e 2 + 0(e 3 ) and hence the motion is sub-diffusive. The time-dependence in 



d = 1 is also known jp5[ , giving the following behaviour for (r ) as a function of dimension 

(r 2 ) oc [log(Dt)} 4 
(r 2 ) oc {Dt) 2 ' z 



(r 2 ) ~ AD R t 




for d = 2 in weak disorder 



(4) 



where, for comparison, the result for pure diffusion is (r 2 ) = 2dDt, and Dr in equation 
effective, measured diffusion constant in the late-time limit. 



is the 



We will show that below two dimensions, a universal form similar to the pure-diffusion d < 2 
result in (||) with n a function of the time t cannot be found. This is due to the changed dynamic 
exponent z which requires a dimensionful amplitude that must be a function of the disorder 
strength or the reaction rate. Though, it is reasonable to study the density as a function of time, 
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or more specifically as a function of the length Dt, it is not the appropriate length scale for the 
reaction-diffusion problem in the presence of disorder. The natural length to use is the disorder- 
averaged diffusion length. By rewriting the density decay as a function of the scale (r 2 ) a fully 
universal relation, similar to (0) can again be found 



" = (^575 with B " 



J_ 210,(12*0-11 

3tt€ 12tt V ; 



for d = 2 - e. 



The effect of Sinai disorder in two dimensions is not strong, and the alteration to the diffusion 
length (H|) is not leading order. Nevertheless, we find the interesting result that a reaction process 
occurring in this disorder has a decay rate with a different leading-order amplitude from the 
MARW, 

24irD R t V 

In fact for weak disorder the reactions occur faster than in the MARW, contrary to the effects 
seen in the case of long-range potential disorder. Nevertheless, when written in terms of the 
diffusion constant of the underlying lattice model (to be described below) it will be shown that 
the disorder strength increases the density's amplitude. These effects come from two competing 
terms: the disorder-renormalisation of the reaction term that increases the rate of reaction and 
the disorder-renormalisation of the propagator that decreases the rate of reaction. 

In the rest of this paper we describe how these results were derived in more detail. In section 
(2) we introduce the model and describe some of the steps taken in the field-theoretic analysis of 
the model. In particular, the relation to existing theories, of diffusion in Sinai disorder |24|| and 
reaction with pure diffusion M || are discussed. The fixed point structure of the renormalised 
parameters is found and a perturbation expansion for the density, valid at early times, is obtained 
for d < 2. In section (3) we obtain a Callan-Symazik (CS) equation for the density as a function of 
time and show that no universal functional form can be found for d < 2. However, by re-expressing 
the density in terms of the disorder-averaged diffusion length a universal form is obtained and the 
amplitude calculated to one-loop order. The behaviour at the upper-critical dimension d c = 2 is 
then examined and the density as function of the disorder strength analysed. Finally, we close in 
section (4) with a discussion of the results obtained. 



2 The model and method 



In this section, we introduce the model to be studied and also describe some of the steps taken 
to achieve its representation in field-theoretic form. The method used is standard and we only 
dwell on details that are different from systems previously studied. After the model has been 
defined it is written in the language of second quantisation, which in turn allows a mapping 
to a path-integral formulation. An average over all possible realisations of the random velocity 
field can be taken at this point, to produce a weighting function (an action) that gives disorder- 
averaged correlation functions. This bare action is then regularised and finally used to calculate a 
perturbation expansion for the early-time, disorder-averaged reactant density. 

The model is defined on an infinite cf-dimensional hypercubic lattice with a lattice spacing of 
unity. Each site i of this lattice contains rii particles where rij can take the values 0, 1, 2 • • •. The 
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quenched disorder in the diffusion rates is modeled by particles hopping independently from a 
lattice site i to a neighbouring site e at a fixed rate Pi—> e - The rates {p} are random and contain 
no long-range correlations. Reactions can occur if there are two or more particles, n, > 2, on a 
lattice site. This happens at a rate Anj(nj — 1), where A is the on-site reaction rate, reducing the 
number of particles on that site by 2. 

The field-theoretic description is obtained by writing a master equation that describes the time- 
dependent flow of probability between microstates. It is convenient to write this equation in the 
language of bosonic operators. Given that the set of occupation numbers {rii} defines a microstate 
of the system, the probability that the system is in such a microstate will be written P({m}). The 
master equation is dt\ip(t)) = —H\tp(t)), where the probability state vector | ?/>(£)) and evolution 
operator 7i are 

\m) = E p (K})IK4) nj i°) 

{ni} j 



n = E 



2 V 2 



(5) 



This algebraic description can now be converted to a field theory by using the coherent-state 
formalism. Observables, like the expected density {rij(t, {p})) at site j at time t for a given 
realisation of the disorder {p} can be written as a path integration with respect to an action iS p 

(nj(t,{ P })) = Jl[ [vfcVtl] fa exp(-S p ). 

i 

The integration is over the complex fields (ft, <ft and the action S p derived from equation (|B|) is 



s P = E 



-4>i{t) + J* dt U\dt4>i + 4 E (P«-e& - Pe^e) ~ A(l - {4>\f)<t>l j " «O0i(O) 



It is convenient to shift the field <f>* by its classical value 4>^ = <ft + 1 and take the continuum limit 
in space. The action S(V) thus obtained is naturally split into four parts Sd + 5^ + Sr + 5 no : 
the diffusive, disorder, reaction and initial conditions. For the moment let us examine the diffusive 
and disorder parts 



S D = j dt j atx (<f)d t (f) - <pV 2 (D(x)0)) 
= - £dtj al d x0V ■ (V(x)cf?j . 



The disorder appears in both the diffusion constant and in a random velocity vector- field V(x). 
However, as can be checked under the RG, the disordered, spatially varying component of D{x) 
is irrelevant in the technical sense, and therefore we consider just the case of a uniform diffusion 
field D{x) = D. Furthermore, though we have used a lattice model in the derivation of the 
continuum field theory, it is not necessary that some of the restrictions of the lattice formulation 
are passed to the continuum theory. In particular, for a walker on a lattice with a bias V there is a 
minimum dispersion D = V/2. In the continuum, there is no reason to impose such a restriction, 
and therefore we treat the magnitude of the diffusion constant and V(x) as fully independent 
quantities. 
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The velocity vector-field V(x) is taken to be a Gaussian random variable with the correlator 
(V a (x)VP(y)) = A6 a>/3 S(x - y), i.e. there are no long-range correlations. An average can be 
performed on the component with respect to this field, to produce an action that gives disorder- 
averaged correlation functions, 

exp (-5 A ) = JVV exp (~ J d d x(V(x)) 2 ^j exp 

This integral can be performed to give the disorder-averaged bare action Sq = Sd + <Sa + <Sr + <S no > 
which when suitably regularised, can be used for calculations. The various components of this 
action are 

rt r _ 



S D = J dt J d d x(j) [d t - D X7 Z ) (6) 

S A = ^ J d d x (J* dt^W^j ■ (J dt^V^j 

S R = 2A f dt f d d x0(f) 2 + A [dt [d d xlf 2j - 2 
Jo J Jo 

S no = —no / d d x / dt(j)5{t). 



o 

Where we have written Dq, Aq and Ao with subscripts to stress that they are bare quantities, but 
no represents the initial density at t = 0. It is seen that the action is a combination of that found 
in p4| for diffusion in Sinai disorder (though with a trivial change of field variables) and that of 
the purely diffusive reaction process ||. Following the notation in previous works, the diagrams 
for the vertices are given in figure (1). The two reaction vertices figure (la), renormalise identically 
and hence we do not introduce a separate reaction parameter for each vertex. It should also be 
noted that the disorder vertex figure (lb), when considered in momentum space, is proportional 
to the scalar product of the momentum flowing through the two outgoing eft fields. Finally, as 
with both the pure reaction and pure disordered-diffusion theories individually, the upper-critical 
dimension for the hybridised theory is d c = 2. 

Because taking the continuum limit in space introduces unphysical divergences, the theory must 
be rendered finite before calculations can proceed. This was achieved by dimensional regularisation 
in d = 2 — e dimensions in the absence of the S no term. In this theory there is no field renormal- 
isation, in fact only the diffusion constant, disorder strength and reaction rate are renormalised. 
The following renormalised diffusion constant D and dimensionless interaction parameters g and 
h are now introduced 

Z D D = D Z g gDn e = Z g X = A Z h hD 2 ^ = Z h A = A . 

It can be shown that the disorder vertex and propagator are not renormalised by the reaction 



vertices, and therefore we can use the results found previously [24] for Zh at the one-loop level and 
Zd at the two-loop level. In fact the only new diagrams that need be considered are the dressings 
of the the reaction strength by the disorder. As can be checked, only one such diagram, see figure 
(1), diverges in d = 2. This contribution is combined with the previously known result for the 
reaction-reaction renormalisation to give the following set of Z factors 

h 2 (g-h) h 

Z D = 1 + jt^t + ■■■ Z g = l+ 1 + ... Zh = l + — + .... 

(4:ir) z e 2ire 4ire 
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(a) reaction vertices (b) disorder vertex 


(d) tree-level density / 
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(e) tree-level propagator # # # 
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(c) renormalisation of a reaction vertex 






(f) density calculation to one loop 




+ _n- + _/- + JiK 

n (0) n (1) ^* % n (1) ^- n (1) * % 

n n g nh n ct 



Figure 1: The various quantities expressed as combinations of the vertices (la) and (lb). The 
requirement that the sum (lc) is finite gives the one- loop contribution to the renormalisation 
constant Z g . Diagram (Id) is the mean-field density given in equation (7). Diagram (le) is the 
propagator dressed by the tree-level reaction interaction and is used in the one-loop contribution 
to the density calculation shown in diagram (If). 



The flow function for the diffusion constant and the beta functions can now be evaluated 



= d\og(D) 
d log(n) 

dg 
d login) 
dh 



2h 2 
(4tt) 2 

a 



+ oih 3 ) 



aiog(/i) 



^ ig-h-2ire) + Oig 6 ,g z h,gh 
-?-(/i-47re)+0(/ l 3 ). 



The fixed point structure for h is of course unchanged by the presence of the reaction vertex, with 
h* = 47re+0(e 2 ) and the dynamic exponent remains z = 2+g* = 2+2e 2 +0(e 3 ). However, the fixed 
point of the renormalised reaction strength is now shifted to the larger value of g* = 6ire + 0(e 2 ) 
in the presence of Sinai disorder (the value with no disorder is g* = 2ire + 0(e 2 )). 



The perturbative density 

With the theory now regularised, observables such as the reactant density n(i) at an early time t 
can be calculated perturbatively 



nit) = J V4>V(j) (p exp i-S - S ct ) 



where the bare action has been rewritten in terms of the renormalised action and the appropriate 
counter term action So = S + S c t. The initial conditions are included and the renormalised 
action is equivalent in form to the bare action written above, but with the replacements Dq — > D, 
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Ao — * hD fi e and Ao — > gDfi* . The early-time density calculation is performed as a loop expansion 
of equation (2) taken to the appropriate order of the small parameter e = 2 — d. An expansion is 
also made in inverse powers of the initial density no with only the leading-order, rig independent 
term retained (the reason for this limit of large initial density will be given in the next section). 

The tree-level or mean-field result shown in figure (Id) is is unchanged from the case of the 
MARW, equation ([[]) 



n 



(o) 



1 



2g^Dt 



(7) 



However, at the one- loop level a new diagram appears, representing two particles interacting with 
the same disordered region at different times, and eventually annihilating. This addition, the third 
shown in figure (If), is equal to 



n 



(i) 



1 



h 



(Dt) d / 2 g 



1 21og(7r)-3 



4ire 



+ 



16vr 



To produce the full, regularised one-loop contribution, this result is combined with the reaction- 



reaction contribution rig , which is also common to ||, and the appropriate counter term n 



(i) 

ct 



n 



(i) 



(Dt) d / 2 



1 | 21og(8vr) 
4ne 16tt 



n 



(i) 

ct 



h-g 
A-Ktg^Dt 



(8) 



The full result can be expanded to finite order in e to produce the full one-loop perturbative result 



n 



(i) 



1 



16vr(L>t) d / 2 



21og(8vr) - 5 - - (21og(7r) - 3) 
9 



+ 0(e). 



(9) 



3 The reactant density 



In this section, the late-time behaviour of the reactant density will be examined. The perturbative 
density, has already been derived and is given by the sum of equations (0) and @. Now a Callan- 
Symanzik (CS) equation will be used to relate this perturbative density to the required late-time, 
non-perturbative density. First, the differential CS equation will be obtained and solved as a 
function of time t and the case of dimensions less than two will then be examined. It will be shown 
that no universal relation giving the density as a function of time can be obtained. However, by 
trading the time dependence in the CS equation for dependence on the disorder-averaged diffusion 
length, a fully universal relation between this length scale and the reaction decay rate will be 
obtained. The behaviour at the upper-critical dimension is then examined. It will be shown 
that the reactant density in two dimensions, written in terms of the measured, late-time diffusion 
constant, decays at a faster rate than the case of pure diffusion-reaction. However, written in 
terms of the diffusion constant of the underlying lattice model, it is shown that as the disorder 
strength increases the reaction rate decreases. 



A Callan-Symanzik equation for the density 

All physical quantities must be independent from the arbitrary \i that was introduced in section 
2 in the definition of the dimensionless interaction strengths g and h. Using this fact, and also 
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dimensional analysis, the following differential equation can be written for the density n 



n 







(10) 



where the interaction dependent flow functions g, (3 g and (3^ are given in the previous section. This 
can now be solved in the standard way by writing equation (|l(]) as a complete differential with 
respect to a scaling variable s, thus 



d 



dlog(s) 



+ d 



n 



0. 



where n = n(s) is the value of the density at some scale. Following the notation used in || we 
denote early-time quantities X(s) (at a scale s) as X and late-time quantities (at a scale s = 1) 
simply as X(l) = X. The flow equations for the system parameters as a function of the scale s 



arc 



dlog(Dt) 



2 + 



Oh 



fa, 



<91og(n ) 



<91og(s) " ' K <91og(s) f " y ' <91og(s) ' <91og(s) 
and the equation relating the density at these two different scales s = 1 and s is 

n(Dt, g, h,no, fi) = s d n(Dt,g,h,no,fi) 



-d. 



(11) 



(12) 



The procedure will now be to insert the perturbative results, (|7|) and (|9|), suitably rewritten with 
g — > g etc, into the RHS of equation (^) and replace all arguments with the s = 1 values using 
the solutions of equations ( |Tl] ) . Before proceeding it should be noted that the solution of the flow 
equation for the initial density implies no = no/s d . As we will be interested in the late-time regime 
s — > 0, this justifies the 1/uq expansion in the perturbative calculation in the previous section. 



Below two dimensions 

The behaviour of the reactant density for d < 2 is now considered. From the flow equation (|I1 
in d = 2 — e it is seen that as s — > the quantities g, h and g approach their fixed-point values 
67re, 47re and 2e 2 respectively. Thus, time varies with the scale s as s z = t/t where z = 2 + 2e 2 at 
the two- loop level [24]. Using these results, the following form for the reactant density below two 
dimensions is found 

"W = w^- < 13 > 

This form appears similar to the d < 2 case in equation (^) but actually lacks all the important 
universal features of that relation. By dimensional analysis, the prefactor C must carry (2 — z)dj z 
units of length and must therefore be a function of A, A or both. No such universal amplitude 
like Ad can be found that relates density to time for the case of reactions in a Sinai velocity field. 
Moreover, given that the dynamic exponent z = 2 + 2e 2 + 0(e 3 ) is itself an approximate quantity, 
there is little point attempting to derive the non-universal prefactor C. However, we will now 
show that, though there is no universality in a relation between the density and the length scale 
Dt in the presence of disorder, a fully universal relation nevertheless exists in terms of the reaction 
density and a different length scale: the disorder-averaged diffusion length (r 2 ) 1 / 2 which is the 
typical distance a single reactant explores. This quantity can be shown to vary as (r 2 ) a t 2 l z , also 
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with a non-universal amplitude. We now rewrite the CS equation ( |i0|) for the density, exchangin 
the dependence in time for a dependence in (r 2 ) 



d\og{r 2 ) 9 dg dh <91og(no) 



d 



n 



0. 



This is solved as above with the scaling s 2 = (r 2 )/(r 2 ). The perturbative density up to the one- 
loop level n^ ) + given in equations (|7|) and (Q) can also be rewritten with the substitution 
2dDt = (r 2 ) which is correct at this order. Combining these results gives the following late-time 
density expressed as a function of (r 2 ) 



1 



n 



■n 



1 

(j,2\jd/2 



(r 2 ) 1 n((r 2 ),5*,/i*,n ,^) 
1 21og(128vr) - 11 



37re 



+ 



12vr 



+ 0(e) 



(14) 



Both n and (r 2 ) are non- universal functions of t, in as much as the disorder strength A enters 
explicitly. However, we have found a universal relation between them, independent of all system 
parameters except the dimension of space. Though the amplitude calculated at this order in an 
e expansion is unlikely to give a good result for one dimension, the scaling relation n oc l/(r 2 ) d / 2 
is exact at all orders in perturbation theory. Therefore, in one dimension it is expected that the 
product n(t){r 2 (t)) 1 / 2 approaches a universal, constant value independent of the disorder strength, 
reaction rate, and initial density. This is in agreement with the result found in [22| for the 
A + A — > O with infinite reaction rate. Given the validity of a factorisation assumption made in 
p2fl for Sinai disorder, the expected exact result for one dimension would have given an amplitude 
1/47T - equivalent to the MARW. 



Two dimensions 

We now consider the case of the reactant density in two dimensions, d = 2. To obtain the 
asymptotic behaviour it will only be necessary to use the tree-level perturbative density 
rewritten in terms of the diffusion length (r 2 ). The relevant equations to be inserted into the 
scaling relation are 

~<o) 2 ~ 67r ( r2 ) 2 

g{r 2 ) log(s) (r 2 ) 

Combining these results yields the following forms for the reactant density in two dimensions 

•*>-£$+° «■*>-') »u=^° (•-)■ <'*> 

where in the second expression the density has been rewritten in terms of time by using the result 



(Q). At this point comparison can be made with the MARW. The result (|2|) and (15) are of the same 
form, but differ from each other by the amplitude: the disorder renormalisation of the reaction 
term has decreased it from l/8n to 1/24-7T. This implies that reactions occur at an increased rate 
in the presence of Sinai disorder. 

It is interesting at this point to consider the behaviour of different lattice models with fixed 
diffusion constant D. By writing a CS equation for the diffusion length the following relation 
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between the effective diffusion constant measured at late time Dr, and the lattice-model parameter 
D can be obtained 



°« " D ( 1 "2^ + ( A2 ))- < 16 > 

The above relation is valid for weak disorder (small A) and implies that the effective diffusion 
constant is reduced from the lattice-model value D by the disorder strength. Hence, if written in 
terms of the long time and length scale behaviour of a lattice model with parameters D and A, 
the reactant density becomes 

In terms of the lattice parameter D it is seen that the reaction rate is still faster than the MARW. 
However, as the disorder strength is increased the reaction rate starts to decrease. Unfortunately, 
the result is valid only for weak disorder and it is not possible to determine from ( [l7[ ) if a point is 
reached where there is a cross-over and the reaction rate becomes less than the MARW. 



4 Discussion 



We have examined the late-time density of reactants in the single-species process A + A — * O in 
the presence of an uncorrelated, quenched random velocity field: so called Sinai disorder. Contrary 
to many existing works on reactions in disorder, the statistics chosen for this velocity field were 
such that there were no long-range correlations. Despite the lack of correlations, it was shown that 
in two dimensions and below disorder changes the density decay on all time scales. The model, 
introduced in section (2), was analysed in the language of field theory and the renormalisation 
group was used to obtain the late-time reactant density. It was shown that the appropriate action 
that generates the disorder-averaged correlation functions is a combination of terms seen in the 
field-theoretic analyses of diffusion in Sinai disorder |24|] and the A+A^O reaction in the absence 
of disorder || . The new interaction diagrams that appear in this hybridised theory were identified 
up to the one-loop level, and the perturbative density was calculated. It was shown that at this 
level a new term appears corresponding to an annihilation of two particles that both interacted 
with the same region of the random field at earlier times. In section (3) the late-time forms for 
the reactant density were derived from the perturbative results, by the use of the appropriate 
Callan-Symanzik equation (|i~0|). 

Below two dimensions the effects of trapping in Sinai disorder are severe and it was found 
that the reaction process becomes sub- diffusion limited, equation Ql3|). It was shown that because 
of the changed dynamic exponent, the relation between the reactant density and time must be 
non-universal, i.e. dependent on the reaction rate or disorder strength. By writing the density as a 
function of time, the pure diffusion length Dt is the implicit length scale. However, in the presence 
of disorder this is an inappropriate scale. Rather, the density should be written as a function of 
the disorder-averaged diffusion- length squared (r 2 ). By exchanging the dependence in t for (r 2 ) at 
the level of the CS equation a universal relation between these two non-universal quantities was 
obtained that is independent of the reaction rate, disorder strength and initial density. 

At the upper-critical dimension d c = 2 the asymptotically-exact form of the density was ob- 
tained as a function of time, equation (|l5|). The random velocity field has the effect of reducing 
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the amplitude of the density decay which implies that for weak disorder the rate of reaction is 
faster than for purely diffusing reactants: an effect coming from the disorder renormalisation of 
the reaction term given in figure (lc). Physically, it represents the process whereby two reactants 
are pushed into the same region of space by the disorder and therefore brought closer together 
than if they were simply diffusing without a bias. Taking the view-point of the lattice model, it 
is appropriate to express the density decay in terms of the model parameter D rather than the 
measured, late-time effective diffusion constant Dr. In this case, it is seen that as the disorder 
strength is increased the reaction rate begins to decrease, equation (|i~7|). This occurs because the 
rate at which the particles explore space is reduced due to the diffusion-constant renormalisation, 
an effect coming from the dressing of the propagator by the disorder vertex. The result obtained 
here is an expansion in the disorder strength A and is therefore only valid only for weak disorder. 
It would be interesting to obtain results for the strong disorder case, perhaps from a numerical 
approach, to see if increasing the disorder further produces a density decay slower than the MARW. 

Finally, we briefly compare the effects of the disorder examined in this paper and the case of 



long-ranged potential disorder |2(| [21|. The effect of potential disorder in two dimensions is more 
drastic because the exponent of the decay is changed, whereas for Sinai disorder the amplitude 
is altered. The relative severity can be understood by the nature of the disordered landscape. 
In a study of diffusion in various forms of disordered landscapes |l3| , it was noted that potential 
disorder produces a landscape with deep trapping wells where, to escape from a trap, any path 
a particle might take involves movement in an unfavourable direction. However, this is not the 
case for Sinai disorder where the landscape (in two dimensions) does not have the morphology of 
potential wells and any pseudo-traps that might exist will tend to have velocity drifts nearby that 
allow for escape. 
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